Hole pockets in the t — J model 
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We present an exact diagonalization study of the electron momentum distribution n{k) small 
clusters of t — J model for different hole concentrations and t/J. Structures in n(fc) which were 
previously interpreted as a 'large' Fermi surface are identified as originating from the well known 
many-body backflow. To obtain reliable information about the true Fermi surface, we focus on 
the regime t<J, where the backflow effect is weak and suppress the formation of a bound state by 
introducing a density repulsion between holes. We find clear signatures of a Fermi surface which 
takes the form of small hole pockets. Comparison of the scaling of n(fe) and that of the quasiparticle 
weight with t/J suggests that these pockets persist also for t>J. 
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I. INTRODUCTION 

The unusual properties of high-temperature supercon- 
ductors have led to great interest in the physics of corre- 
lated electrons near a Mott-Hubbard metal-to-insulator 
transition. Thereby a particularly intriguing problem is 
the volume of the Fermi surface (FS) for the slightly less 
than half-filled band: should one model the doped insu- 
lator by a dilute gas of quasiparticles corresponding to 
the doped holes (this would imply that the volume of the 
FS is proportional to the hole concentration) or do all 
electrons take part in the formation of the Fermi surface, 
so that its volume is identical to that of free electrons? 
It is the purpose of this paper to present evidence that 
for finite clusters of t—J model the first picture is the 
correct one: the FS as deduced from the momentum dis- 
tribution takes the form of small hole pockets. The t—J 
model reads: 

H--t E {£lc^,. + H.c.) + jJ2[S.-S,-^]. 
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The Si are the electronic spin operators, cl ^ = —cl^{l — 
Ui^-a-) and the sum over < i,j > stands for a sum- 
mation over all pairs of nearest neighbors. Various au- 
thors 0-0] have computed the momentum distribution 
naik)={cl. ^Ck,a) for the two-hole ground state of small 
clusters of this model (corresponding to a nominal hole 
concentration of ~10%) and found it roughly consis- 
tent with a free-electron picture: n{k) is maximum at 
fc=(0,0), minimum at fc=(7r,7r). It has become cus- 
tomary B to cite this as evidence that already at such 
fairly low hole concentrations the t—J model has a free 
electron-like ('large') FS. It is straightforward to see, 
however, that this shape of n{k) is simply the conse- 
quence of elementary sum-rules and has no significance 



for the actual topology of the FS 0]. We have there- 
fore performed a systematic study of the n{k) for various 
doping levels and t/J. 



II. SINGLE HOLE CASE 

As compared to the uniform value of 1/2 for the 
half-filled case, the introduction of only a single hole 
changes n{k) in a rather complex way. Fig. n^ shows 
n{k) for the single hole-ground states with momentum 
^0 = (7i'/2, 7r/2) in the 16-site cluster and momentum 
fco — (27r/3, 0) in the 18-site cluster. The k dependence 
of n(A;) is roughly consistent with free electrons, i.e. n{k) 
is large near (0,0) and small near (7r,7r). This structure, 
which simply ensures negative kinetic energy W, is less 
pronounced the smaller t/ J . The second characteristic 
feature are 'dips' at fep for the minority spin (i.e. the 'hole 
spin') and at fco + (7'',7r), for both spin directions. These 
dips are more pronounced for smaller t/ J . The question 
arises which of these features should be associated with 
the FS, i.e. do we have a 'large' FS already for a single 
hole or is there a 'hole pocket' at fco? We note that the 
magnitude of the discontinuity in n{k) has to be equal 
to the weight of the quasiparticle peak in the single par- 
ticle spectral function, Zh. Since Zh has a pronounced 
Q (and therefore characteristic) dependence on t/J, a 
potential FS discontinuity must have the same charac- 
teristic dependence on t/J. Then, the 'depth' of the dip 
at fcp can be estimated by comparing with a symmetry 
equivalent fc-point i.e. for fco = (7r/2,7r/2) we consider 
^d^p = ni(-7r/2,7r/2) - n^(7r/2, 7r/2), for fco = ( W3, 0) 
we study A^ip = nj^(0,27r/3) — n4(27r/3, 0). In Fig. g 
these differences are compared to Zh (obtained from the 
single particle spectral function for momentum transfer 
fco at half-filling) for various values of t/ J . Obviously, 
Arfij, = Zh over the entire range of t/J, so that that the 
dip clearly originates from the Fermi level crossing of the 



quasiparticle band, i.e. we have a 'hole pocket' at feg- On 
the other hand, differences An(fc) across the 'large' FS 
always show the opposite behaviour under a variation of 
t/ J as Zh, indicating that these drops in n(fc) are unre- 
lated to any FS crossing. This suggests to associate this 
structure in n[k) with the well-known 'backflow' for in- 
teracting Fermi systems [pi . Such a strong backflow effect 
is by no means surprising if we consider the change in the 
single particle Greens function upon removing one elec- 
tron with momentum fco near the FS: naively one might 
expect that the only effect be the shift of the quasipar- 
ticle peak at fco from the photoemission to the inverse 
photoemission spectrum. If this were true, however, the 
integrated photoemission weight (which equals the to- 
tal number of electrons) had decreased only by Z/j ^ 1. 
Hence the bulk of spectral weight shift from photoemis- 
sion to inverse photoemission must occur for momenta 
fc 7^ fco, i.e. the strong backflow. What is remarkable 
is the wide spread of the backflow in fc-space. This im- 
plies that for each individual fc the change of n{k) due 
to removal of an electron at kp is ~ 1/-/V, with N the 
system size. For a small finite hole concentration 5 it 
seems reasonable that the backflow contributions from 
the individual holes are additive, so that the total back- 
flow contribution would scale with 5. 
What remains to be explained are the 'satellite dips' at 
fco + (tt, tt). The most natural explanation are antiferro- 
magnetic spin correlations. To see this, let us consider 
the case i <C J, where the state l/\/2cka,i\'^o) (with |$o) 
the half-filled ground state) to good approximation is an 
eigenstate. For this state 



{nk.x) 



9 V k,ko' 



ls,k 



fco), 



14 1 

S'(fc-fco) + — , 

2 3 ^ "^ iV' 



(1) 



with S{q) the static spin structure factor of |$o)- Since 
the latter is peaked sharply at q — (tt, tt) we have a nat- 
ural explanation for the 'satellite dips', and it seams rea- 
sonable to adopt this explanation also for larger values 
oit/J. 

Summarizing the results obtained so far, we may say that 
the introduction of a single hole changes n{k) in a rather 
complex way: there is a dip at the momentum of the 
hole, which originates from the Fermi level crossing of 
the quasiparticle band and thus represents the 'Fermi 
surface'. The dip is superimposed over a smooth free- 
electron-like variation, the familiar many-body backflow. 
As a consequence of the small quasiparticle weight, this 
backflow is very pronounced in the t — J model. For later 
reference we note that the backflow contribution to n(fc) 
to good approximation is a function of \kx\ + \ky\ only 
(this is also confirmed by investigating n{k) for other 
fco). Finally, the strong antiferromagnetic spin correla- 
tions produce dips also at ko + Q. 

This shape of n{k) can be easily understood by recall- 
ing M that the elementary excitations near the Fermi 
energy are spin bags, where the hole is dressed by anti- 



ferromagnetic spin fluctuations; a simple calculation in 
terms of the string picture ||^ reproduces the numerical 
results quantitatively. 



III. TWO HOLE CASE 

We proceed to the ground state with two holes. Var- 
ious authors ||l|-01 have found that the free electron-like 
variation of riffc) observed already for a single hole be- 
comes more pronounced for this doping level, and based 
on the criterion n{k) > 1/2 M the 'Luttinger Fermi sur- 
faces' in Fig. would be assigned. However, by the same 
arguments as for a single hole, these Luttinger Fermi sur- 
face are ruled out: Fig. ^ compares the t/J dependence 
of differences An{k) across the respective Luttinger FS 
to that of the quasiparticle weight in the spectral func- 
tion for the two-hole ground state. Zh decreases sharply, 
the An(fc) increase monotonically with t/J. The drop in 
n(fc) upon crossing the large FS thus is obviously unre- 
lated to any true Fermi level crossing. Instead, compar- 
ison with Fig. g shows that the t/J dependence of the 
An(fe) is very similar to the backflow contribution for a 
single hole. More precisley, if we assume that the back- 
flow for the two holes is simply additive, we expect for 
the 'large FS' differences in the two-hole ground state: 
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where An^ (fc) are the corresponding differences in the 
single hole ground state. Fig. || compares the 'large FS' 
An(fe) with the estimates obtained from (||) by using the 

Aria (fc) shown in Fig. |2[ Both the magnitude and the 
t/J scaling are predicted very well by (0), which clearly 
suggests to associate the 'large FS' with the backflow 
contribution. 

The question then is: what is the true FS and how can 
we make it visible in n(fc)? Fig. |l| shows that for a single 
hole the true Fermi surface (i.e. the hole pocket ) is 
most clearly visible for large J/t. This is simply related 
to the fact that the quasiparticle weight is large in this 
parameter region (see Fig. ||). Since the i/J-scaling of 
Zh is essentially the same at half-filling and in the two 
hole ground state, (see Figs, g and Q) we thus may expect 
to see the clearest FS signatures for large J/t also in the 
two-hole ground state. 

For more than one hole, however, we face an additional 
problem: the strong interaction between the holes, which 
manifests itself e.g. in a sizeable negative binding energy 
MEOl. An interacting state of two 'quasiparticles' reads 
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Thus, whereas for a single hole we could fix the location 
of the pocket simply by choosing the total momentum, 
in the two-hole ground state the holes will be distributed 
over different momenta with probability ~ |A(fc)p and 



one may not hope to observe any FS signature unless 
A(A;) is well localized in fc-space, i.e. A(A;) ~ Sk,ko with 
the quasiparticle ground state kg. This in turn neces- 
sitates that the interaction energy be smaller than dif- 
ferences in single particle energy between neighboring k- 
points, i.e. weak interaction and sufficiently strong dis- 
persion. With this in mind we add a density interaction 
term Hv=V ^^^ ,y niUj, to the Hamiltonian; adjusting 
the parameter V, one may hope to reach a situation, 
where Hy to a certain degree 'cancels' the intrinsic at- 
tractive interaction of the holes. In addition, we include a 
small next-nearest neighbor hopping term in the Hamil- 
tonian, so as to lift the unfavourable (near) degeneracy 
of the quasiparticle dispersion along the surface of the 
magnetic Brillouin zone; we fix the value of the respective 
hopping integral to be t' = — O.lt. For the 16-site cluster 
this term has the additional advantage that it breaks the 
spurious additional symmetry due to the mapping to a 
2^ hypercube and selects a unique two-hole groundstate 
with momentum (0,0). 

Let us stress the following: due to the addition of these 
terms we are strictly speaking no longer considering the 
original t—J model. It seems quite plausible, however, 
that if a FS exists at all, its volume should be changed 
neither by changing the kinetic energy (i'-term) nor by 
introducing an additional interaction (X^-term). 
To demonstrate the adjustment of V, Fig. shows 
the variation of the hole density correlation function 
g(i2)=(l/2)2,((l - nK,)(l - hr^+r)) in the two-hole 
ground states of the 16 and 20 site cluster with V (due 
to a subtle but understandable pathology in its geometry, 
analogous results cannot be obtained for the 18-site clus- 
ter, see Appendix). Its essentially identical behaviour in 
both clusters clearly signals a change of the net interac- 
tion between the holes from attraction to repulsion. For 
intermediate values oiV, on the other hand, g{R,) is quite 
homogeneous indicating that the single particle delocal- 
ization energy of the holes dominates over their interac- 
tion. Given this plus the large Zh we should therefore be 
in an optimal position for observing the FS. 
Figs. ^ and g show the single particle spectral function 
A(fe,w) for J/t=2 and the momentum distribution for 
J/t=l and J/t=2 (and the respective 'optimal' V). In 
the spectral function, the chemical potential Ep is lo- 
cated near the top but within a group of pronounced 
peaks, well separated from another such group in the 
inverse photoemission spectrum. There are pronounced 
peaks both immediately above and below Ep which com- 
prise the bulk of spectral weight for the respective mo- 
menta. Corresponding to the well defined 'quasiparticle 
peaks' in the spectral function, n(fe) exhibits a sharp vari- 
ation: hole pockets at (7r,0) and (0,7r). They are super- 
imposed over the familiar backflow contribution, which 
again has the generic free electron like form so as to en- 
sure negative kinetic energy. Fig. | also gives the val- 
ues of the quasiparticle weight for the 'Fermi momenta'. 
For (27r/5, Stt/S) in the 20-site cluster the 'quasiparti- 
cle peak' in the phtoemission spectrum (PES) actually 



consists of two peaks with approximately equal weight; 
we consider these as a single 'broadened' peak, so that 
the two weights should be added. This is supported by 
the good agreement with Zh at (7r/2, 7r/2) in the 16-site 
cluster (where no splitting occurs) and the reasonable 
agreement with the Zh deduced from the inverse pho- 
toemission spectrum (IPES) at (tt, 0) (where no splitting 
occurs either). The 'depth' of the pockets approximately 
equals Zh and both quantities consistently decrease with 
decreasing J/t. 

A somewhat astonishing feature of these results is the lo- 
cation of the pockets at (tt, 0) rather than at (7r/2, 7r/2). 
This can be traced back to the point group symmetry of 
the ground two-hole ground state: when the symmetry 
of the half- filled ground state is Ai (or s), that of the 
two-hole ground state is Bi (or d^2_y2) and vice versa 
(the former situation is realized in the 16 and 18-site 
cluster, the latter in the 20-site cluster). Addition of two 
holes thus always is equivalent to adding an object with 
dj.2_j^2-symmetry, which implies that the pair wave func- 
tion A(fc) in (0) should have this symmetry as well. This 
in turn implies A(fc) = for fc along (1, 1), so that occu- 
pation of (tt, 0) is favoured. 

Fig. pi shows the 'FS discontinuities' in the 4x4 cluster 
under a variation of the repulsion strength V . They show 
maxima when the density correlation function is most ho- 
mogeneous, precisely as one would expect for Fermions 
with a variable interaction strength. In the spectral func- 
tion the reduction of the discontinuities as y — > man- 
ifests itself by the reduction in intensity of the big IPES 
peak at (tt, 0) and the appearance of small low energy 
IPES peaks at the momenta next to (7r,0): the pockets 
are 'washed out'. We note that the An{k) across the 
'large' FS remains unaffected, another indication that it 
is unrelated to low-energy physics. 
A possible explanation for the hole pocket FS would be 
spin-density-wave-type broken symmetry: although the 
ground states under consideration are spin singlets, this 
might be realized if the fluctuations of the staggered mag- 
netization Ms were slow as compared to the hole mo- 
tion, so that the holes move under the influence of an 
'adiabatically varying' staggered field. A possible crite- 
rion for this situation would be Ttr-i^AF'^'^''^, where Ttr 
is the time it takes for a hole to transverse the cluster 
and lliaf is the frequency of fiuctuations of Ms- We 
estimate (for J/t = 2) the group velocity of the holes 
from the dispersion of the 'quasiparticle peak' in the 
PES spectrum and, using the energies indicated by ar- 
rows in Fig. u\ for the 20-site cluster and the peaks at 
{tt/2, 0) and {tt/2, 7r/2) for the 16-site cluster, we find 
Ttr^27r/0M{2TT/0.2t) for the 20 (16)-site cluster. Typical 
frequencies for fluctuations of Ms can be obtained from 
its correlation function, which, up to a constant, equals 
the dynamical spin susceptibility for momentum transfer 
(tt, tt); a rigorous lower bound on loaf thus can be ob- 
tained by subtracting the ground state energy from the 
energy of the lowest state with total momentum (tt, tt) 
and the same point group symmetry as the ground state. 



This gives w^F>0.9t(1.2i) for the 20 (16) cluster, i.e. 
Ttr ■ ujaf>'^t^- 'Almost static' Neel order thus can be 
ruled out as origin of the small FS, even for this fairly 
large value of J/t. 

As an additional check we have introduced exchange 
terms J' between 2"'* and 3'"'' nearest neighbors to re- 
duce the spin correlations and again optimized the re- 
pulsion to enable 'free' hole motion. Ground state prop- 
erties of this (highly artificial) model are summarized in 
Fig. nffl: the momentum distribution, hole density cor- 
relation function and spin correlation function S'di?!) = 
exp{iQ ■ R){Si ■ Si+R) (with Q — (tt, tt)). The density 
correlation function is homogeneous (no charge ordering), 
the spin correlations decay rapidly (no long range anti- 
ferromagnetic or spiral ordering) but still there are un- 
ambiguous hole pockets in n(fe). The only possible con- 
clusion is that it is only the large Zh which makes the 
pockets visible in the large J region, and not the onset 
of any kind of ordering. 

While the hole pockets can be made clearly visible for 
large J, the situation is more involved for t>J. In this 
parameter region the small overlap between 'quasiparti- 
cle' and 'bare hole' (as manifested by the small Zh) makes 
the y-term (which couples only to the bare hole) in- 
creasingly inefficient in enforcing a noninteracting state: 
rather than separating from each other, the two holes 
remain bound on second-nearest neighbors up to fairly 
large values of V and the 'crossover' from attraction 
to repulsion, which gave an unambigous prescription for 
choosing 1/, cannot be obtained any more. We thus aban- 
don both the V and t' terms and adopt a more indirect 
way of reasoning. 

In the single hole case, we found that the pocket was su- 
perimposed over the smooth backflow contribution. We 
assume that the situation for 2 holes is similar, only with 
the additional complication that the pockets are now 
'washed out' due to the interaction between holes (see 
Fig. 1^). Therefore we expect that n{k) can be written 



n{k)^nback{k) + Zh-\Hk)'[' 



(4) 



with the pair wave function A(A;) introduced in (0). As 
discussed above, the point group symmetry of the two- 
hole ground state necessitates that A(fc) has dx^-y2 sym- 
metry, so that the pockets are located at (7r,0). Then, 
since the symmetry of the ground state is unchanged by 
adding either the t' or V term, we conclude that this 
should also hold true in the absence of these terms. 
Since n{,ack{f^) to good approximation is a function of 
\kx\ + \ky\ only (see Fig. P, this contribution can be 
eliminated by forming the difference of two momenta 
with (almost) equal |fcx| + |fcy|- Next, if we choose one of 
these momenta along (or near) the (1, 1) direction, where 
the (i^2_j,2 -symmetry requires that A(fc) vanishes (or is 
small), and the other at (or near) (tt, 0) we should obtain 



An{k)^Zh-\A{Tr,0)\' 



(5) 



so that, in contrast to the 'large FS' differences indicated 
in Fig. y this difference should scale with Z^- To check 
this prediction, the t/J dependence of various such dif- 
ferences is shown in Fig. |ll| and obviously, they are to 
excellent approximation proportional to Z^ over a wide 
range of t/J. The scaling of n{k) with t/J thus is com- 
pletely consistent with the assumptions 

a) that there are washed out hole pockets at (7r,0), 

b) that these are superimposed over the smooth backflow 
contribution, which is the sum of the backflows for the 
two individual holes (see Fig. 0). 



IV. COMPARISON WITH OTHER NUMERICAL 
STUDIES AND EXPERIMENT 

While the hole pockets are very clearly visible for large 
J, the evidence for their existence in the physical regime 
is of a more indirect character. A comparison with 
other numerical calculations and experiments on high- 
temperature superconductors is therefore necessary. As 
far as numerical studies on small clusters are concerned, 
hole pockets and/or rigid band behaviour upon doping 
are consistently suggested by most of the available nu- 
merical calculations. For the t ~ J model, Poilblanc and 
Dagotto |ll| studied the A{k, u;) for single hole states and 
concluded that the two-hole ground state in the 4x4 clus- 
ter shows hole pockets at (tt, 0), in agreement with the 
present result. On the other hand, Stephan and Horsch 
studied n{k) and A{k, uj) for the two-hole ground state 
and concluded that there is neither rigid band behaviour 
nor hole pockets. However, these authors based their con- 
clusions solely on the qualitative inspection of a rather 
limited data set, which is largely irrelevant M for de- 
ciding the FS topology. In addition to the inconsistent 
scaling behaviour found above (Fig. 0) , numerical calcu- 
lation of A{k,uj) for the 20-site cluster [Q rules out the 
Luttinger FS postulated by Stephan and Horsch. 
Castillo and Balseiro [|2| computed the Hall constant 
and found its sign near half-filling to be consistent with 
a hole-like FS, i.e. with hole pockets. Gooding et al. M] 
studied the doping dependence of the spin correlation 
function in clusters with special geometry and also found 
indications of rigid-band behaviour. Finally, a system- 
atic study of the doping dependence of the single par- 
ticle spectral function [|l4| shows rigid-band behaviour, 
i.e. holes are filled into the quasiparticle band present at 
half- filling (which naturally implies hole pockets). 
The situation is quite similar for the Hubbard model. 
While the generic Q free-electron like shape of n{k) 
found in earlier Monte-Carlo studies ||l^ was initially 
considered as evidence against hole pockets, more careful 
and systematic analysis [ |l6| showed that hole pockets are 
in fact remarkably consistent with the numerical data, 
their nonobservation in the earlier studies being simply 
the consequence of thermal smearing. It seems fair to say 
that the available numerical results for small clusters of 
both Hubbard and t~J models, when interpreted with 



care, are all consistent with rigid band behaviour and/or 
hole pockets. 

Let us next discuss experimental results on high- 
temperature superconductors assuming that the hole 
pockets found in the cluster studies persist in the real 
systems. The volume of the FS associated with the 
Cu — O plane-derived bands in these materials presents 
a well-known puzzle: early photoemission experiments 
|_8| show bands, which disperse towards the Fermi en- 
ergy and vanish at points in /c-space which are roughly 
located on the free electron FS corresponding to elec- 
tron density 1 — S, where S is the hole concentration; 
on the other hand transport properties can be modelled 
well ||l^,^ by assuming a FS with a volume ~ S. In 
a Fermi liquid, the apparently contradicting quantities 
actually fall into distinct classes: photoemission spectra 
depend on Zh, transport properties do not. Hence, if one 
wants to resolve the discrepancy entirely within a Fermi 
liquid-like picture, the simplest way would be to assume 
a 'small' FS and explain the photoemission results by a 
systematic variation of Z^, along the band which forms 
the FS, similar to the 'shadow band' picture pl[. A 
trivial argument for such a strong fe-dependence of the 
quasiparticle weight is, that a distribution of PES weight 
in the Brillouin zone (and hence a n(fc)) that resembles 
the nointeracting FS, always optimizes the expectation 
value of the kinetic energy. Therefore it is favourable 
if those parts of the band structure, which lie inside the 
free-electron FS have large spectral weight, and the parts 
outside small weight. Then, it seems that in a recent pho- 
toemission study by Aebi et al. H^ , structures which are 
very consistent with such a shadow band scenario have 
indeed been observed. Moreover, another key feature of 
the dispersion relation for a single hole, namely the ex- 
tended flat region near (tt, 0) p3,pj| has also been found 
as an universal feature of high temperature superconduc- 
tors i2a,M Adopting a rigid band/hole pocket scenario 
thus would explain many experiments in a very simple 
and natural way, which moreover is remarkably consis- 
tent with the existing numerical data as a whole. 



V. CONCLUSION 

We have discussed the problems in directly determin- 
ing the FS from n{k) in the t—J model: small quasipar- 
ticle weight, a pronounced 'backflow' effect and strong 
interaction between the doped holes. We have then ex- 
amined the single particle spectral function and n(k) in a 
situation where these problems were largely avoided and 
found signatures of a FS which takes the form of small 
hole-pockets. Analysis of the scaling of n{k) with t/ J 
suggested that these pockets also persist in the regime 
t > J . Available numerical data all support this picture, 
and we have outlined a possible scenario to reconcile ex- 
periments on high-temperature superconductors. 
The assumption of a small Fermi surface implies that 
the phase of some given basis state is determined by a 



Slater determinant of rank N — N^, {N^ being the num- 
ber of electrons) rather than Ng (as it would be e.g. in a 
Gutzwiller projected Fermi sea). Moreover, it would not 
be the positions of the electrons which enter this Slater 
determinant, but those of the hole-like quasiparticles, so 
that we have a very different nature of long range phase 
coherence. The Fermi surface in an interacting system, 
being a 'remnant' of the noninteracting one, is obviously 
a consequence of the requirements to have minimum ki- 
netic energy and to satisfy the Pauli principle. On the 
other hand, close to half-filling most of the electrons are 
immobile, so that the gain in kinetic energy from creating 
the long range phase coherence between electrons (which 
is responsible for the singularity of n{k)) may not be 
very large. On the other hand, the vacancies are almost 
unconditionally mobile, so that phase coherence between 
holes may be more favourable. 
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VI. APPENDIX 

The 18-site cluster has a pathological geometry, which 
does not allow for an unbound state of two particles with 
d2,2_j^2-symmetry. The ultimate reason is that the prim- 
itive lattice translations in this cluster are (3, 3) and 
(3,-3). Writing an interacting two-hole state in real 
space, we have 
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If we choose R — (2, 1), rotate counterclockwise by 7r/2, 
reflect by the a;-axis and add (3, 3) we recover the orig- 
inal vector. A state with d^-i^yi symmetry picks up a 
factor of (—1) during these operations, hence 0(2, 1) = 0. 
Analogous reasoning shows that (^(3, 0) = 0, so that all 
large distances between particles are 'symmetry forbid- 
den' (possible distance in this cluster are (1,0), (1,1), 
(2,0) (2,1) and (3,0)). The 'unbinding transition' for 
this cluster thus can occur only via level crossing and 
this is indeed the case: when V is switched on in the 
d^2_y2 ground state, the holes stay close to each other 
even for fairly large values of V. Instead, at V ^ 3t (for 
J/t = 2) a level crossing occurs, and a new ground state 
with momentum (27r/3,0) is stabilized. 



FIG. 1. Momentum distribution for the single hole ground 
state with Sz = 1/2 (i.e. with a '|-hole') of the 4x4 cluster 
with t/J = 4 (top) of the 4x4 cluster with t/J = 1 (middle) 
18-site cluster with t/J = 2 (bottom). The upper values refer 
to the majority spin, the lower values to the minority spin, 
the ground state momentum fco is marked by a black box and 
fco + (tt, tt) by a dotted box. 



FIG. 10. Hole density correlation function g{R) (dark cir- 
cles), spin correlation function S{R) (light squares) and n{k) 
(inset) for the ground state with spin frustration. Parameter 
values are J = 2, J' = 0.75, t' = —0.1. There is a density 
repulsion of strength 1.7t between holes on 1°*, 2"'' and 3'^'' 
nearest neighbors. 



FIG. 2. Comparison of the i/J-dependence of Zh at 
half-filling (dark squares) and various differences An(fe) in 
the single hole ground states with Sz ~ 1/2. Shown is the 
'depth' of the pockets (light circles) and differences across the 
'large FS' (up and down triangles). 



FIG. 11. Comparison of the scaling of Zh (dark squares) 
and selected difference 

An(fe) (light squares) with t/J: 1.5-(n(7r/2, 7r/2) -n(7r,0)) 
(16-site), 2.8-(n(7r/3, 7r/3) - n(27r/3,0)) (18-site) and 
3.6-(n(7r/5,37r/5) -n(7r,0)) (20-site). 



FIG. 3. Allowed momenta and 'Luttinger Fermi surfaces' 
for various clusters. The 'Fermi momenta' are denoted by kp- 



FIG. 4. Comparison of the t/J-dependence of Zh (dark 
squares, obtained from the photoemission spectra at the kp 
indicated in Fig. H) and the differences An(fe) between pairs 
of momenta connected by dashed lines in Fig. t| (light squares 
and circles). 



FIG. 5. Comparison of the t/J-dependence of the An(fc) 
shown in Fig. g (light squares) and the estimates obtained 
from (pi) by using the single hole An(fc) in Fig. EI (dark 
squares) . 



FIG. 6. V-dependence of the hole density correlation func- 
tion g{R) in the ground state of t—t' — J—V clusters with 2 
holes (J/t=2). 



FIG. 7. Single particle spectral function for the t—t' — J—V 
model with 2 holes and J/t=2 for the 20-site with V/t=2.5 
(a) as well as the 16-site cluster with V/t=2.4 (b). Delta func- 
tions have been replaced by Lorentzians of width 0.05t. The 
frequency region lu<Ef {lu>Ef) corresponds to photoemis- 
sion (inverse photoemission). 



FIG. 8. Momentum distribution in the two-hole ground 
state of the cluster t~t'~J~V models with J/t—2 and 
V/t=2.5 {V/t=2A) for the 20-site (16-site) cluster (a) and for 
J/t=l and V'/i=3.0 (V/t=2.0) for the 20-site (16-site) cluster 
(b). For the 'Fermi momenta' the quasiparticle weight Zh is 
given in brackets. 



FIG. 9. V dependence of selected An(fc) in the 4x4 cluster, 
parameters are like in Fig. H. 
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